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We find the constraints from WMAP and SDSS data on the fraction of cosmological fluctuations 
sourced by local cosmic strings using a Markov Chain Monte Carlo (MCMC) analysis. In addition to 
varying the usual 6 cosmological parameters and the string tension (/r), we also varied the amount 
of small-scale structure on the strings. Our results indicate that cosmic strings can account for 
up to 7 (14)% of the total power of the microwave anisotropy at 68 (95)% confidence level. The 
corresponding bound on the string mass per unit length, within our string model, is Gp < 1.8(2.7) x 
10“^ at 68 (95)% c.l., where this constraint has been altered from what appears below following the 
correction of errors in our cosmic string code outlined in a recent erratum, astro-ph/0604141 0. We 
also calculate the B-type polarization spectra sourced by cosmic strings and discuss the prospects 
of their detection. 

PACS numbers: 98.80.Cq 


NOTE: The results of the analysis presented in this 
paper are somewhat modified in light of the discovery of 
errors in the computer code we use in this paper. These 
errors, and their effects on the results, are discussed in 
a new erratum to this paper, astro-ph/0604141 Q. The 
text of this paper has not been altered. Please read it 
together with our erratum to learn what is still valid and 
what has changed. 


I. INTRODUCTION 

Cosmic strings were proposed as potential seeds for 
structure formation in the early 1980’s Hi and exten¬ 
sively studied in this context until the late 1990’s (see 0 
for a review). The interest was fueled, in part, by the 
fact that the CMB temperature anisotropy discovered 
by COBE 0 was of the same order of magnitude (10“®) 
as the density fluctuations that would be produced by 
cosmic strings formed around GUT epoch. However, by 
1997 it was already apparent that strings could not ex¬ 
plain the distribution of large-scale structure |^. Finally, 
when Boomerang Q and Maxima revealed the exis¬ 
tence of acoustic peaks in the CMB angular spectrum, 
strings were excluded as a viable alternative to inflation 
as a model for seeding the formation of structure. The 
shape and location of the peaks, currently measured to 
a high precision by WMAP strongly indicate that 
structures grew out of initial fluctuations with adiabatic 
initial conditions and with a scale-invariant spectrum, 
just as prescribed by the inflationary paradigm M . In 
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contrast, strings predict nearly featureless CMB spectra^ 
and, when normalized to COBE or WMAP, do not pro¬ 
duce enough clustering on large scales. 

While strings could not have seeded all of the struc¬ 
ture in the Universe, they could have created a subdom¬ 
inant yet non-negligible fraction of the primordial cos¬ 
mological fluctuations 0. This 

idea has recently received renewed attention with the 
realization that cosmic strings are produced in a wide 
class of string theor y m odels of the inflationary epoch 
mmilllllllllllllllll. In these models, in¬ 
flation can arise during the collisions of branes that co¬ 
alesce to form, ultimately, the brane on which we live 
mils mm. Brane inflation predicts adiabatic tem¬ 
perature and dark matter fluctuations capable of repro¬ 
ducing all currently available observations. In addition, 
the collision at the conclusion of brane inflation can pro¬ 
duce a network of local cosmic strings mm , whose 
effects on cosmological observables range from n egligibl e 
to substantial, depending on the specific scenario |2,‘lil24l| . 
It has also been shown that strings could form at an ob- 
servationally acceptable level at the end of the D-term 
inflation in SUSY GUT models [T^. As the precision of 
cosmological observations increases, one might hope to 
distinguish among the numerous presently-viable models 
of inflation by studying and constraining the properties 
of the cosmic strings they predict. 

The properties of the strings produced at the end of 
brane inflation are similar to those of local strings. They 


^ This is true for all “conventional” cosmic strings networks. Tech¬ 
nically, it is possible to construct active sources of fluctuations 
leading to CMB temperature spectra in agreement with observa¬ 
tions, if, for example, one allows the time- and space-c orrelati ons 
of strings to violate the bounds imposed by causality [iniii 
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intercommute and break off loops like networks of usual 
strings, implying that they will settle into a scaling so¬ 
lution. However, because of their higher-dimensional na¬ 
ture, they are expected to intercommute at a reduced 
rate and thus approach a scaling solution at a higher 
string density. There are also other, qualitatively new 
aspects of cosmic superstrings, the most important of 
which is that string theory models generically predict the 
existence of at least two fundamentally different kinds of 
cosmic strings. These are the so-called fundamental, or 
T’-strings, and another kind of string formed by wrapping 
of all but one of the dimensions of a higher-dimensional 
Z^-brane around compact dimensions, or iZ-strings. The 
relative tensions of these string types are determined 
by the superstring coupling, gs- Because these different 
string types interact through binding, rather than inter¬ 
commutation, they can form higher-tension (p, q) states 
composed of p F-strings and q F-strings yj] ■ A first 
model of these networks has been developed recently, and 
shows that these binding reactions allow network scaling 
1^. Because of the complexity introduced by this kind 
of new physics, the amount of small-scale structure on 
strings - the string “wiggliness” - can vary depending 
on the particular brane inflation scenario. To make our 
constraints applicable to a wide range of models, we have 
included the string wiggliness as one of the parameters in 
the MCMC simulation. The effect of a higher string den¬ 
sity can be roughly modeled by appropriately adjusting 
the normalization of the string-generated spectra. This 
then translates into a simple adjustment of the bound on 
Gp. 

The aim of this paper is to constrain the properties 
of cosmic strings by using the power spectrum data from 
the WMAP and SDSS experiments. There are other ways 
to constrain cosmic strings - some of them promising to 
produce much tighter bounds than those that will ever 
be possible with power spectrum data. We give a brief 
review of other methods in the summary section, Sec. nvi 

The rest of this paper is organized as follows. In Sec. QH 
we give a detailed account of the model and the methods 
used. We show the results in Sec. cni and conclude with 
a discussion in Sec. CYI 


II. THE MODEL AND ANALYSIS 

The fluctuations resulting from brane inflation are ex¬ 
pected to be an incoherent superposition of contributions 
from adiabatic perturbations initiated by curvature fluc¬ 
tuations and active perturbations induced by the decay¬ 
ing cosmic string network. The resulting CMB angular 
spectra can be written as a sum of the adiabatic and 
string contributions: 

Ci = C^’^ + Cf\ ( 1 ) 

Analogous expressions hold for matter density spectra. 
We restrict our study to a flat FRW universe and vary 



10 100 1000 
1 


FIG. 1: The CMB TT and TE spectra (solid lines) sourced 
by cosmic strings with wiggliness parameter Ur = 1.9, as well 
as the adiabatic spectra for the same cosmological parameters 
(dashed lines) and WMAP’s first year data. The string spec¬ 
tra are normalized so that the total TT power is the same for 
the two lines, which corresponds to B = 1. 


the following cosmological parameters: the Hubble con¬ 
stant h, the matter density the baryon density 

and the reionization optical depth r. In addition, 
we vary the galaxy bias factor Ff,, the amplitude and 
the spectral index Ug of the primordial scalar power spec¬ 
trum, as well as the string mass per unit length, G/i, 
and the string wiggliness parameter (to be defined in 
Sec. Ill All . 

We used a suite of different codes to produce and ana¬ 
lyze the spectra. The model we employed for the cosmic 
string-generated perturbations is described in the sub¬ 
section below. The string CMB and matter spectra were 
calculated using a modification of CMBFAST 

(see Figs. ^ and |2 for representative string induced 
spectra). We first evaluated and stored the string spec¬ 
tra on a grid of parameters. During our calculations, the 
spectra were obtained by interpolation on the grid. The 
adiabatic matter spectra were also stored on a grid after 
having been evaluated using a publicly available version 
of CMBFAST. For the adiabatic CMB spectra, we used 
the package CMB Warp [^ . 

To compare the theoretical linear matter power spec¬ 
trum P^{k) generated by our code with the SDSS results, 
we first applied the halo-fitting procedure of Smith et al. 

to obtain the non-linear spectra, P^^{k). This pro¬ 
cedure is only valid for the adiabatic contribution to the 
P{k) spectrum, so using it for the cosmic string contri¬ 
bution introduces some inaccuracy into our model. How¬ 
ever, as we can see in Fig. 12 the string power spectrum 
is considerably weaker than the adiabatic power spec- 
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FIG. 2: The string-generated matter power spectrum (solid 
line) for the same parameters as in Fig. 0 i.e., for B = 1. 
The dashed line represents the linear spectrum from adiabatic 
perturbations at 2 = 0. 

trum on all but the smallest length scales even in models 
where strings account for all the large-scale microwave 
anisotropy. Thus, since only a small portion of the rele¬ 
vant values of k are affected by the halo-fitting procedure, 
and since the effect is small on the scales we consider, we 
expect any inaccuracy introduced by our approximate 
treatment of nonlinearity to be unimportant. Next, we 
assumed that the galaxy bias factor and the correction 
due to peculiar velocities on the scales of our interest 
{k < 0.2 h/Mpc) are approximately redshift- and scale- 
independent and can be combined into a single constant 
factor, Ff,, multiplying the matter spectrum: 

P{k) = FbP^^ik). (2) 

We then fit P{k) to the SDSS power spectrum data (for 
k < 0.2 h/Mpc) with Fb as a free parameter. 

The likelihood of the spectra for a given set of param¬ 
eters was obtained using the publicly-available likelihood 
codes produced by the WMAP ^3 SDSS ^3 teams. 

A. The Cosmic String Model 

Unlike the adiabatic inflationary perturbations, which 
are set as initial conditions in the equations of motion 
of linear perturbations, cosmic strings act as a contin¬ 
ual source of fluctuations as the universe evolves. The 
density and temperature fluctuations created in the im¬ 
mediate locality of a cosmic string are highly non-linear, 
e.g. the density contrast 6p/p is significantly larger than 
unity in the wake formed behind a moving string. How¬ 
ever, the effect of strings on cosmological scales is that of 


a small perturbation to the evolution of the average cos¬ 
mic energy-momentum tensor. Hence, for the purpose 
of computing the CMB and LSS spectra, the metric and 
density perturbations caused by strings are described by 
the system of linearized Einstein-Boltzmann equations 
with strings acting as active sources. 

Evaluation of the CMB and LSS spectra sourced by 
strings requires knowing their energy-momentum tensor 
(or its unequal time correlation functions ) for the en¬ 
tire dynamical range of the calculation, which is approxi¬ 
mately four orders of magnitude in scale factor. Realistic 
simulations of cosmic string networks have so far been 
limited either in their dynamical range or their res¬ 
olution ^3- Hence, until full scale simulations become 
available, one is forced to resort to approximate methods 
to model the string sources. 

Numerical simulations [dHl Id^ show that during 
the radiation and matter dominated eras the string net¬ 
work evolves according to a scaling solution, which on 
sufficiently large scales can be described by two length 
scales. The first scale, ^{t), is the coherence length of 
strings, i. e. the distance beyond which directions along 
the string are uncorrelated. The second scale, L{t), is the 
average inter-string separation. Scaling implies that both 
length scales grow in proportion to the horizon. Simula¬ 
tions indicate tha t £(f) ^ t, while L{t) = yt, with y » 0.8 
in the matter era |^|^. The one-scale model |^|d3, in 
which the two length scales are assumed to be equal, has 
been quite successful in describing the general properties 
of cosmic string networks inferred from numerical sim¬ 
ulations. These simulations assume that cosmic strings 
reconnect on every intersection. It is of interest to us, 
however, to also consider the case when the reconnec¬ 
tion probability is less than one. If strings can move and 
interact in extra dimensions then, while appearing to in¬ 
tersect in our three dimensions, they may actually miss 
each other. Hence, the effective intercommutation rate 
of these strings will generally be less than unity. As a 
consequence, one would expect more strings per horizon 
in these theories MM Because of the straightening 
of wiggles on sub-horizon scales due to the expansion of 
the universe, one would still expect ^(t) ^ t. However, 
the string density would increase, therefore reducing the 
inter-string distance. Hence, smaller inter-commutation 
probabilities imply smaller y. 

In addition, numerical simulations show that long 
strings possess a great deal of small-scale structure in 
the form of kinks and wiggles on scales much smaller 
than the horizon. To an observer who cannot resolve this 
structure, the string will appear to be smooth, but with 
a larger effective mass per unit length p, and a smaller 
effective tension T. An unperturbed string (with p = T) 
exerts no gravitational force on nearby particles. In con¬ 
trast, a wiggly string with p > T attracts particles like a 
massive rod. The effective equation of state of a wiggly 
strings is pT = p"^ [ 5 ( 1 . 1 ^ . Depending on the brane infla¬ 
tion model, the presence of extra dimensions could mean 
that even more small scale structure would be present on 
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the strings [^. 

To calculate the sources of perturbations we use an 
updated version of the cosmic string model first intro¬ 
duced by Albrecht et al.@ and further developed in refs. 
[^ 1 ^, where the wiggly nature of strings was taken into 
account. In this model, the string network is represented 
by a collection of uncorrelated straight string segments 
produced at some early epoch, moving with uncorrelated, 
random velocities. At every subsequent epoch, a cer¬ 
tain fraction of the number of segments decays in a way 
that maintains network scaling. The length of each seg¬ 
ment at any time is taken to be equal to the correlation 
length of the network. This length and the root-mean- 
square (r.m.s.) velocity of segments are computed from 
the velocity-dependent one-scale model of Martins and 
Shellard |^. The positions of segments are drawn from 
a uniform distribution in space, and their orientations are 
chosen from a uniform distribution on a two-sphere. This 
model is a rather crude approximation of a realistic string 
network. However, there are good reasons to believe that 
its predictions for the CMB and LSS spectra are close to 
what one would obtain by using full-scale simulations. Its 
parameters have been calibrated to produce source cor¬ 
relation functions in agreement with those in [H^ . where 
comparison to a full simulation was possible. Also, the 
shape of the spectra obtained using this model are in 
good agreement with results of other groups , 

who used different methods that are also approximate. 
Finally, on the cosmological scales probed by the CMB 
measurements, the fine details of the string evolution do 
not play a major role. It is the large-scale properties - 
such as the scaling distance, the equation of state (wig- 
gliness), the r.m.s. velocity, and how all these charac¬ 
teristics evolve through the radiation-matter equality - 
that determine the shape of the string-induced spectra. 
All of these effects are accounted for in our model and 
can, in principle, be adjusted to match any specific cos¬ 
mic string model. We choose to work with this model 
because one can easily calculate the sources for different 
cosmological parameters and because it allows us to in¬ 
clude the effect of the wiggliness , which could be one 
of the distinguishing features of strings produced in the¬ 
ories with extra dimensions 1^. The other main effect 
of the presence of extra dimensions, the increased string 
density, can be approximately factored in by multiplying 
the spectra by Ns ~ 7 “^. 

For technical details of the model, the reader is re¬ 
ferred to |3f)| and references therein. The wiggly nature 
of strings is accounted for by modifying the string energy- 
momentum tensor so that it corresponds to the wiggly 
string equation of state: 

fl = apL , T = , (3) 

where a is a parameter describing the wiggliness, p and 
T are the mass per unit length and the string tension of 
the wiggly string, and fi is the tension (or, equivalently, 
the mass per unit length) of the smooth string. In addi¬ 
tion to modifying the equation of state, the presence of 


small-scale structure slows strings down on large scales. 
We account for this by dividing the root mean squared 
string velocity by the parameter a. The wiggliness of 
the strings remains approximately constant during the 
radiation and matter eras, but changes its value during 
the transition between the two. We take the radiation era 
value, ar, to be a free parameter that we vary, and set the 
matter era value to be am = (l+ccr)/ 2 , with a smooth in¬ 
terpolation between the two values (as prescribed in [s^). 
For conventional strings, this roughly agrees with results 
of numerical simulations [iM which show a decrease 
from ar ^ 1.8 — 1.9 in the radiation era to am ~ 1.4 — 1.6 
in the matter era. 

The way the wiggliness modifies the shape of the spec¬ 
tra can be understood qualitatively, once one examines 
the main physical processes that contribute to CMB 
anisotropy on various scales. The spectrum can be ap¬ 
proximately divided into two main components. One is a 
roughly flat (scale-invariant) component that arises from 
the combined Kaiser-Stebbins effect of many strings 
in a scaling string network; this component is produced 
after last scattering. The second component is the bump 
peaked at £ ^ 450, which is primarily determined by the 
state of density and velocity perturbations at the time of 
last scattering. CMB photons experience Doppler shifts 
when they (last-)scatter off the velocity flows created by 
the string wakes that exist at this epoch |^. In ad¬ 
dition, all the density fluctuations that existed at the 
time of last scatterin g ar e imprinted on the CMB via 
the Sachs-Wolfe effect ; these density fluctuations are 
caused by the CDM wakes created by strings during the 
time between radiation-matter equality and last scatter¬ 
ing. If we change the string substructure by introducing 
string wiggliness, there is a two-fold effect on the flat part 
of the spectrum that arises from the Kaiser-Stebbins ef¬ 
fect. The first effect is simple: by nature, heavier strings 
have larger deficit angles, and thus produce larger CMB 
discontinuities. However, the size of the discontinuities 
is also proportional to string velocity and, since wiggly 
strings move more slowly than straight ones, the net in¬ 
crease in the size of the discontinuities is mitigated. The 
wakes behind a wiggly string are, nonetheless, always 
more prominent. This is because of wiggly strings hav¬ 
ing a non-zero Newtonian potential that attracts matter 
particles in the string’s vicinity. The combined effect of 
wiggliness on the flat component of the spectrum is an 
overall increase in its amplitude, but a relatively weak 
one. The peak, on the other hand, is directly enhanced, 
since the wakes at the last scattering are more prominent 
for wiggly strings. The net result is an overall increase 
in the amplitude of the spectra coupled with an enhance¬ 
ment of the peak at £ ^ 450. The characteristic scale of 
the fluctuations at decoupling is set by the size of the typ¬ 
ical wake, or, in turn, by the coherence length ^(£) of the 
string network at that time. Independent of the string 
intercommutation probability, the coherence length ^(£) 
is always expected to be roughly equal to t. This is be¬ 
cause of the straightening of wiggles on subhorizon scales 
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caused by the expansion of the universe. Therefore, the 
qualitative features of the spectrum are quite indepen¬ 
dent of the details of the model. Hence, any observa¬ 
tional constraint that we obtain based only on the shape 
of the spectrum (such as the constraint on the parameter 
/ defined below) can be viewed as less mo del-dependent 
than our other results. 

It is well known that properties and possible obser¬ 
vational signatures of global and local strings can be 
very different. Global strings predict almost no power 
on small angular scales for CMB temperature anisotropy 
[ 4 ^ . while local strings - as we have argued above - 
produce a quite significant broad peak at I ^ 450 in a 
spatially flat universe |3^ |5^ |^; this can 

be seen in Fig. Q] Global strings also induce a signifi¬ 
cantly larger vector component of metric perturbations 
1^. Consequently, their prediction for the strength of 
the H-type polarization is generally higher than that 
of local strings. We will only concern ourselves with local 
strings. For the most recent constraints on global strings 
the reader is referred to |6^ and [^. 

Rather than working directly with the string parame¬ 
ter G/i, we introduce a parameter B, defined as 



where po is the tension that one obtains by setting the 
total power in string induced CMB anisotropy to be equal 
to the total power observed by WMAP. That is, we set 

= E , po) = , (5) 

I 

where we take 0 ^°^ = 1-9 and po is a fixed set of the 
remaining cosmological parameters that have taken to 
correspond to the the best fit ACDM model of Tegmark 
et al 1^. The value of Gpo that we obtain with this 
prescription is 2 x 10“®. Speaking very loosely, B can 
be said to measure the fraction of the anisotropy due to 
strings. Note, however, that this meaning is modified 
if, for instance, the strings have a different amount of 
small-scale structure (ar ^ 1.9) or if they have reduced 
intercommutation probabilities. B is only really useful 
as an intermediate parameter. Our main results are the 
constraints on Gp that we obtain from B and the fraction 
/ of the total CMB anisotropy due to strings (Eq. 0). 

B. Markov Chain Monte Carlo 

Because of the large size of our parameter space (nine 
parameters in total), we have used the Markov Chain 
Monte Carlo (MCMC) method for exploring the likeli¬ 
hood surface and for generating marginalized posterior 
distribution functions for the model parameters. We em¬ 
ployed the MCMC algorithm described in the Appendix 

of ii. 



TABLE I: Prior constraints on the parameters 


Parameter 

B = 0 

B > 0 

/ 

— 

< 0.068 (68%), < 0.14 (95%) 

B 

— 

< 0.029 (68%), < 0.062 (95%) 

Ol.j' 

— 

< 2.3(68%), < 3.6 (95%) 

As 

n q'7+0-08 

U.o/ _0.i6 

0.851^3 

Us 

1 

-*-•'^-0.04 

1.0 ±0.026 

Q.Bh^ 

0.024 ± 0.001 

U.UZO_Q 0016 


0.15 ±0.01 

n 1 

U.10_o.oi 

h 

0.69 ± 0.03 

0.71 ±0.034 

T 

0.155 ±0.057 

0.143 ±0.054 

Fb 

J-."±o_Q 22 

1 

-*-•^'-0.18 


TABLE II: The best fit results 


We ran eight separate chains initialized at randomly 
generated initial positions within our prior range. The 
priors, given in TableHl were chosen with the expectation 
of the string contribution being subdominant and the 
best fit parameters being close to their WMAP best fit 
values 1^^. Since we expected values of B near zero to 
be preferred, we also allowed B to range slightly below 
zero so as not to restrict artificially the ability of the 
chain’s random steps to explore near B = 0; values below 
zero were discarded when the data were analyzed. One 
advantage to our use of multiple chains rather than a 
single, long chain was that we were able to verify directly 
that there was adequate mixing in each chain, since each 
successfully forgot its starting location and located the 
same maximum likelihood region of the parameter space. 


^ The restriction on wiggliness, ar < 4, was to save computing 
time. 
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FIG. 3: The one-dimensional projected PDFs for the 9 pa¬ 
rameters varied by our Markov Chain Monte Carlo code; note 
that lob = , u)m = Q,Mh? ■ The solid line represents 

the PDFs for models where cosmic strings are included; the 
dashed line represents the PDFs for models with B = 0, i.e. 
without cosmic strings. Each curve has been rescaled such 
that its area is unity. For each PDF the lightly shaded re¬ 
gions are excluded at the 68% confidence region; the dark 
regions are excluded at the 95% confidence region. 


III. RESULTS 
A. The fraction in strings and 

To test our MCMC code we first ran a chain with the 
string contribution set to zero (i? = 0). The results are 
shown in Table im Our results are consistent with those 
found in a similar analysis of the same data by members 
of the SDSS team [^. 

Each of our eight chains allowed variation in all nine 
parameters. The results are summarized in Fig. 0 and 
Tableim In Fig. Owe plot the marginalized 1-D posterior 
distribution functions for the nine parameters we allowed 
to vary in our analysis. The solid lines represent the 
PDFs for these parameters with cosmic strings included; 
the dashed lines show the results without cosmic strings 
(5 = 0). We have shaded the regions excluded at 68% 
(light) and 95 % (dark) confidence. The peaks of each 
of these PDFs and the one-sigma error bars are given 
in Table HU for the parameter r, where the results lack 
a clear peak, we have taken the midpoint of the 68 % 
confidence region as our “peak” value. 

The resulting cosmology with cosmic strings included 
is very close to the cosmology without cosmic strings. 
This verihes our hypothesis: a subdominant admixture of 
cosmic strings into the cosmological perturbation spectra 
gives a minor modification to the resulting cosmological 


FIG. 4: The PDF for the Cosmic String weighting parameter, 
B = (G/J./2 X 10 -®)^ 


parameters as determined in such a model. This justi¬ 
fies the approximation made in 1 ^ , where we fixed 
all of the cosmological parameters except for Us to their 
WMAP best-fit values. The only noticeable shift caused 
by including cosmic strings comes in the peak likelihood 
in the PDF for and even this shift is small. The 

only potentially worrisome aspect of these PDFs is the 
fact that the PDF for r appears to be running into the 
upper value of our prior range. 

Let us examine the PDF for the cosmic string weight¬ 
ing parameter, 5, given in Fig. The light and dark 
shaded regions again represent the 68 and 95% confi¬ 
dence intervals. Although the par ameter name, B, is re¬ 
tained from our previous work |69| , the fact that we have 
changed our string model has changed this parameter’s 
meaning. Thus, Fig. should not be directly compared 
with our earlier results. Indeed, we have added a new 
parameter for describing the strings, a,., the string wig- 
gliness. In our previous code, we had effectively fixed 
ar = 1. Since strings with ar > 1 have a larger effective 
mass per unit length than ar = 1 strings, the cosmologi¬ 
cal perturbations caused by > 1 strings are generally 
larger. Note that ar and B were not degenerate parame¬ 
ters - the effect of a higher ar-, for a particular spectrum, 
was quite independent of 5; see a (rough) contour plot 
in Fig. O This is because varying the wiggliness modifies 
the shape of the spectra in addition to their amplitude, 
as we described briefly in smin. 

To quantify the total string contribution to the CMB 
anisotropy for a given set of parameters s we define the 
fraction / as 

TCS 

= JCS + jad ( 6 ) 
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FIG. 5: A contour plot of a separate set of MCMC results 
in the two-dimensional space defined by the parameters B 
and a. The contours are somewhat crude, but suggest that a 
is poorly constrained. At this level, there is no evidence for 
degeneracy between the parameters. 


where 


and 


jcs ^ , 


47r 




where are the temperature (TT) correlation spec¬ 

tra. We then compute the PDF for this new parameter, 
which we show in Fig. El along with the 68 and 95% 
probability regions. We find that about 7%(14%) of the 
CMB power can be sourced by strings at 68%(95%) con¬ 
fidence. Note that the dashed line, corresponding to the 
three parameter analysis of |69| , shows a smaller allowed 
string contribution. 

Figs. 0 and El are the principal results of this paper. 
Limits on B alone give limits on the string tension itself. 
Using the results in Fig. El we find a cosmic string weight 
of B < 0.029(0.06) allowed at the 68 (95)% confidence 
level. This corresponds to Gfj, 3.4(5.0) x 10“^. The 
peak of the PDF for B lies aX B = 0.01, or G/r ~ 2x 10“^. 
These limits are relevant to searches for direct detections 
of cosmic strings, as the magnitudes of gravity wave and 
lensing events caused by cosmic strings depend directly 
on their tension. 

The above bounds on G/i are model-dependent. In or¬ 
der to match the total observed CMB power, our string 
model requires G/tq ~ 2 x 10“®. This value is on the up¬ 
per end of other estimates in the literature, which means 


that our upper bound on G/r can be treated as a con¬ 
servative upper bound. Our G/io is consistent with the 
COBE normalized values in (G/r ~ 1.7x 10“®), in 
(G/r = n^±0.5)xl0-®),inl7| (G^ = (1.7±0.7)x 10“®), 
and in (G^ 2 x 10“®). It also agrees with (for 

similar model parameters). Sig nificantly different esti¬ 
mates of G/io were found in [^, where the COBE nor¬ 
malized value of G/i was ~ l.Ox 10“® (for their parameter 
= 1/3), and similarly in |^. The latest estimates 
of Landriau and Shellard using realistic simulations 
of cosmic strings (reliable up to ^ ^ 20) give the 
COBE normalized value of G/i = (0.74 ± 0.2) x 10“® 
for a ACDM cosmolo gy, which is consistent with results 
of a similar study in |74l where the value obtained was 
G/i = (1.05 ± 0.3) X 10“®. Note that most of the results 
obtained prior to 1999 assumed a CDM dominated cos¬ 
mology - switching to ACDM leads to a ~ 10% increase 
in COBE normalized value of G/i |4^ . 

Our bound on G/i would also be altered if the strings 
intercommute at a rate less than unity, as is expected in 
many string theory models of cosmic strings. The effect 
of reduced intercommutation would be to lower the upper 
limit on G/i. 

Our bound on the fraction of CMB power in strings, 
/, depends only on the shapes of the string-sourced CMB 
and LSS spectra. These shapes, as discussed in im are 
largely independent of the details of the string model. 
The bound on / can be used to derive an approximate 
bound on G/i, given one’s favorite value of /iQ. For exam¬ 
ple, if one accepts the values in Hi III, i.e. G/io ~ 10®, 
one obtains G/i' ~ /^^^G/io ^ 2(3) x 10“^ at 68 (95)% 
confidence level. 

It is also worth recognizing that / can serve as a 
measure of the goodness-of-fit of the paradigmatic in¬ 
flationary scenario in comparison with a physically moti¬ 
vated model; isocurvature is another example of a model 
used in such a manner (e.g. [Z^). Our results from 

cosmic strings, serving from this viewpoint merely as a 
self-consistent foil to the standard model, show that as 
much as 14 % of the CMB TT-correlation power could 
be sourced by a radically different spectrum without de¬ 
stroying the close agreement of the resulting spectrum 
with the anisotropy data. Loosely speaking, we can con¬ 
clude from this bound that there is a cumulative am¬ 
biguity in the uniqueness of the adiabatic Ci spectrum 
shape, as determined from the WMAP data, of around 
10%. As more CMB data become available in the future, 
repetition of this analysis might be worthwhile, if only 
to discover whether the intrinsic shape of the adiabatic 
spectrum is more uniquely picked out by the more-exact 
future data sets. 


B. The B-polarization spectrum 

In Fig. 0 we plot the B-mode polarization spectra in 
the case of smooth strings {ttr = 1) and wiggly strings 
Or = 1.9 predicted by our string model for the case when 









Relative String Contribution, f 


FIG. 6: The PDF for the combination parameter /, which 
quantifies the fractional contribution of cosmic strings to the 
total spectrum. The solid line shows / from our full 

analysis, with the 68 % (light) and 95 % (dark) confidence 
regions shaded. The dotted line is the result for / from our 
previous three parameter analysis. 


the total contribution of strings to the CMB anisotropy is 
10%. That is, for each of the curves, the value of G/i was 
adjusted separately to correspond to / = 0.1. For com¬ 
parison, we also plot the B-mode spectra from a purely 
adiabatic cosmology. The light dotted line represents the 
B-mode polarization arising from gravitational lensing of 
E-mode polarization. The light dash-dotted line repre¬ 
sents the B-mode arising from gravitational waves. Any 
direct detection of cosmic-string generated B-modes will 
rely upon accurate predictions of the spectrum of B-mode 
polarization arising from lensing of E-mode polarization. 
Fortunately, several observational groups believe that it 
will be possible to “clean” as much as 90 % of the appar¬ 
ent E to B lensed power through accurate E-mode obser¬ 
vations. With gravity-wave generated B-modes peaking 
at a much lower £, any excess power in observed B-mode 
spectra at £ ^ 1000 could be a telling sign of cosmic string 
activity. Two planned experiments, QUIET and QUaD 
[t^ , expect to be able to measure such high-£ polarization 
with great precision. It is also worth noting that the am¬ 
plitude of gravity-wave generated B-mode polarization 
is intimately tied to the as-yet undetermined scalar-to- 
tensor ratio, r. We have used r = 0.1 in Fig. |7| which 
is usually regarded as an optimistically high value. Infla¬ 
tionary models frequentl y p roduce orders-of-magnitude 
lower estimates for r (in |^, for instance, investigations 
predict 10“® ^ r < 10“^ for KKLMMT-motivated brane 
inflation). For r <C 0.1, cosmic strings could be the dom¬ 
inant source of B-mode polarization for low values of I 
(10 < i < 1000), but with a spectrum that is recogniz- 



FIG. 7: The B-type polarization spectra for ar = 1.9 (dashed 
lined) and a^. = 1 (solid line) both corresponding / = 0.1. 
The light dotted line is the B-mode expected from gravi¬ 
tational lensing of adiabatic E-mode polarization; the light 
dash-dotted line is the B-mode expected from gravity wave- 
sourced polarization, for tensor-to-scalar ratio r = 0.1. 

ably distinct from the polarization generate d by gravity 
waves. The proposed CLOVER experiment and its 
space-based successors plan to focus their measurements 
on this region of Gspace. 


IV. SUMMARY 

We found two types of constraints on cosmic string 
networks through our analysis of the WMAP and SDSS 
data. One is on the value of G/r, which is sensitive to 
our normalization convention and to the string intercom¬ 
mutation rate, which we take to be unity in setting our 
bound. The other is on the fraction of total CMB power 
due to cosmic strings, /. The constraint on / depends 
chiefly on the general shape of the string induced spec¬ 
trum. The shape of the spectrum is quite generic - a 
plateau on large scales due to superposition of the Kaiser- 
Stebbins effect of many strings (from a scaling network), 
followed by a broad peak on small scales. The peak is 
caused by the Doppler and Sachs-Wolfe effects produced 
by velocity and density perturbations caused by strings 
during the epoch of recombination. The rough position 
of the peak is set by the size of a typical wake at last 
scattering, while the wake size is set by the coherence (or 
curvature) length ^{t) at that time. We can call these fea¬ 
tures generic because they are agreed upon by all groups 
that have studied string-induced CMB spectra. The fine 
details of the shape, e.g. the sharpness of the peak, de¬ 
pend on many factors. Among the most important of 
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those extra factors is the amount of small-scale struc¬ 
ture on strings (the wiggliness), which we included as a 
parameter in our calculation. 

In addition to quantifying the allowed fraction of the 
cosmic string contribution, our parameter / can also be 
interpreted as a measure of the goodness-of-fit of the fidu¬ 
cial adiabatic CMB spectrum model. The fact that the 
data permit approximately 10% of the TT-correlation 
power to arise from a very different competitor model 
gives a useful measure of how uniquely the data pick out 
the shape of the adiabatic spectrum. 

Other recent constraints on local strings that also used 
WMAP data include ref. where the narrowness of 
the first peak was used to constraint the size of the inco¬ 
herent string contribution, and refs. where the 

constraint was based on the expected non-Gaussian sig¬ 
natures induced by strings. Ref. suggests an interest¬ 
ing way to obtain a rough bound the string contribution 
to the CMB spectrum at ~ 220 from constraints on 
the width of the main peak. Our method has the advan¬ 
tage of including the information on all scales (not just 
the main peak) and the ability to account for changes in 
in cosmological parameters and the shape of the string 
spectra (by varying the wiggliness). 

Our limit on the string tension, G/x < 5 x 10“^, does 
not contradict a variety of recent claims of observational 
evidence for the existence of cosmic strings. The most 
prominent have been possible examples of cosmic string 
lensing. In the first case [sj, a pair of nearly identi¬ 
cal images were found, the best explanation of which 
appears to be a direct lensing by a cosmic string with 
tension G/x ~ 4 x 10“^. Another possible cosmic string 
lensing observation is the appearance of short time- 
scale variations in the brightness of the well-known grav¬ 
itational lens system, Q0957 -I- 561, which could be ex¬ 
plained by a passing cosmic string loop with a tension in 
the range 10“® < Gfj. < 6 x 10“^. The inferred time scale 
of these variations is so short, ~ 100 days, that very small 
scale string structure would be required, making the 
claim somewhat problematic (for more on these events, 
see [13 )■ It is also possible that strings are responsi¬ 
ble for the early reionization suggested by WMAP [^ . 
since cosmic string wakes can bring about star formation 
much earlier than the standard Press-Schechter sce¬ 
nario conservative analysis of string-mediated early 
reionization suggests a bound of Gfj. < 5 x 10“”^ [sj- In 
the arena of gravitational wave observation, current pul¬ 
sar timing bounds are still marginally consistent with 
G/x ~ 10~^ [S^, while analysis of gravitational wave 
bursts from string loops suggests that G/x ^ 10“”^ cos¬ 
mic strings will be readily observable by both LIGO and 
LISA 3 

It is often thought that a key observational test for cos¬ 
mic strings would be whether non-Gaussianity is found in 
the primordial perturbations seen by such experiments as 
WMAP. This is a natural assumption, since each string 
acting alone would perturb the CMB in a highly non- 
Gaussian manner. However, the central limit theorem 


tells us that the superposition of perturbations produced 
by many strings must be Gaussian. Therefore, one ex¬ 
pects to see string-related non-Gaussian features only on 
scales that are sufficiently small not to have been crossed 
by more than a few strings during the entire period of 
time during which strings have produced their effects. 
It is not difficult to get a conservative estimate of this 
scale. The dominant contribution to the anisotropy on 
small (sub-degree) scales comes from the Doppler and 
Sachs-Wolfe effects produced at the last scattering |^. 
A natural length scale to start with is the angular size 
of the horizon at recombination, which corresponds to 
£ ~ 220/-\/3. Numerical simulations show that 

the typical distance between strings during matter dom¬ 
ination is L ~ 0.8t. It is somewhat smaller in the radia¬ 
tion era, and can be much smaller if the intercommuta¬ 
tion probability was less than unity, as may be the case 
in string theory models. A conservative estimate of the 
number of strings per horizon at any time in the mat¬ 
ter era is ^ (47rt^/3)/L^ ^ 10. When projected onto 
the last scattering surface, about half of this number of 
strings would contribute. Hence, the GMB anisotropy in 
a patch corresponding to £ ~ 220/-\/3 would be a super¬ 
position of the effects of about 5 strings, and to isolate 
the effect of one string, one would have to go to scales of 
order: £ ~ 220 x 5/-\/3 ~ 600. In doing this rough esti¬ 
mate we have ignored the density perturbations created 
by wakes between the radiation-matter equality and last 
scattering, which contribute a non-negligible fraction of 
the power near the main peak |90| . This contribution 
would tend to make the scale at which non-Gaussianity 
appears even smaller, since the wakes started to from 
from the onset of matter domination. Thus one will 
likely need a resolution of at least £ ~ 1000 to have 
any hope of seeing non-Gaussianity from strings. Fur¬ 
thermore, in the above argument we did not account for 
the fact that strings can produce only 10% of the to¬ 
tal anisotropy. This makes the detection of their non- 
Gaussian signatures even more difficult. The possibility 
of low string intercommutation rates (high string den¬ 
sity), in addition to making strings more Gaussian (via 
the central limit theorem), also strengthens the bound on 
their tension, hence further complicating the detection 
of their non-Gaussian properties. The analysis in refs. 
[7^ 113 assumed £ ~ 200 as the scale for the onset of 
non-Gaussianity in the CMB caused by strings. A more 
realistic scale, as we have argued above, is likely to be 
an order of magnitude smaller, so the detection of string 
sourced non-Gaussianity appears to be beyond the reach 
of WMAP, and quite possibly even Planck. The exist¬ 
ing constraints on string-sourced CMB non-Gaussianity, 
such as those obtained in refs. appear to reflect 

the limitation that, given the variance of a CMB map on 
a certain scale, cr^, one naturally has difficulty resolving 
any detailed features on those scales that have amplitudes 
comparable to In light of this, it is not surprising that 
the upper bound on Gfi obtained in [731131 ~ 10“®) 

roughly corresponds to the variance of the WMAP CMB 
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map on sub-degree scales. 

We find that cosmic strings with tensions of Gfi = 
5 X 10“^ are still allowed by the data from the WMAP 
and SDSS experiments, and that strings can account for 
as much as 14% of the the temperature fluctuations in 
a cosmic microwave background radiation dominated by 
adiabatic fluctuations without any significant changes in 
the underlying cosmology. Indeed, the excess power at 
very small angular scales seen in CMB observations may 
already suggest the presence of cosmic string-generated 
perturbations, which may dominate on such scales |9l| . 
Strings with allowed tensions are produced in brane in¬ 
flation models, implying that such models are still vi¬ 
able and that the strings produced by them may be 
observable, giving us hope of an observational window 
on string theory. One promising signature of cosmic 
strings with these tensions in the early universe would 
be their creation of observable B-mode polarization in 
the CMB with spectra distinct both from those created 
by E-mode lensing and by gravity waves from primor¬ 


dial tensor modes. If primordial tensor modes are weak, 
cosmic strings could be the principal new physics seen 
by B-mode observations; successful observation of such 
B-modes would in turn be our first direct observational 
probe into the physics of string theory. 
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